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The structure of homogeneous turbulent shear flow is studied using data generated by Direct 
Numerical Simulations (DNS) and a linear analysis for both compressible and incompressible cases. 
At large values of the mean shear rate, the Rapid Distortion Theory (RDT) limit is approached. 
Analytical solutions are found for the inviscid compressible RDT equations at long times. The 
RDT equations are also solved numerically for both inviscid and viscous cases. The RDT solutions, 
confirmed by the DNS results, show that the even order transverse derivative moments of the dilata- 
tional and solenoidal velocity fields are anisotropic, with the dilatational motions more anisotropic 
than their solenoidal counterparts. The results obtained for the incompressible case are similar to 
those obtained for the solenoidal motions in the compressible case. The DNS results also indicate 
an increase in the anisotropy of the even order transverse derivative moments with the order of the 
moment, in agreement with the RDT predictions. Although the anisotropy decreases with Reynolds 
number, it is likely that for higher even order moments it will persist at large values of the Reynolds 
number, in contrast with the postulate of local isotropy. The RDT solutions also predict that the 
normalized odd order transverse derivative moments of the solenoidal velocity for the compressible 
case and of the velocity for the incompressible case should approach a constant different than zero 
at large times. This prediction is supported by the DNS data. For higher odd order normalized mo- 
ments, the RDT analysis suggests that the anisotropy may persist at large values of the Reynolds 
number, in agreement with the experimental data. The amplification of the dilatational kinetic 
energy in the direction of the mean shear and the anisotropy of the dilatational dissipation tensor 
found in the DNS results are also consistent with the RDT analysis. 

PACS numbers: 47.27.Ak,47.27.Eq,47.27.Gs,47.40.-x 

I. INTRODUCTION 

Most turbulent flows are anisotropic at large scales. Homogeneous shear flow represents one of the simplest 
anisotropic flows and its study can reveal important aspects of the structure and production of the turbulent fluctua- 
tions. The high Reynolds number experiment of Shen and Warhaft 1 indicates that the higher odd order moments of 
the velocity derivatives may not be consistent with the postulate of local isotropy, which requires that the normalized 
odd order transverse derivative moments approach zero at large Reynolds numbers. This finding has important con- 
sequences, since this postulate has been central for turbulence theories and models^ Although DNS at high Reynolds 
numbers are not feasible yet, persistent anisotropy in the skewness of velocity derivatives has been observed earlier in 
the numerical results of Pumir^ also confirmed by Schumacher^ 

The assumption of local isotropy also requires relations between the even order transverse derivative moments, which 

should satisfy ^ (§fj) y = ^ (if^) y ' an< ^ similar relations for 112 and U3. There is experimental evidence that the 

second order transverse derivative moments are not isotropic, especially when the mean rate of strain is significant 
Nevertheless, the anisotropy of the higher even order transverse moments of the velocity derivatives has not been 
examined. In this study we point out a lack of isotropy of the even order transverse derivative moments in homogeneous 
shear flow for both incompressible and compressible cases with increased anisotropy for higher order moments. We 
also show that the small scale dilatational motions are more anisotropic than their solenoidal counterparts, so that 
the anisotropy increases for the compressible case. Moreover, the anisotropy of the even order transverse derivative 
moments found in the DNS results is predicted by the RDT solutions. As the Reynolds number increases, the 
anisotropy among the transverse derivative moments is expected to decrease. Since it is shown that the anisotropy 
increases with the order of the moment, some anisotropy may persist for higher order moments at large values of the 
Reynolds number. In addition, it is shown that the linear analysis predicts that the normalized odd order transverse 
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derivative moments of u\ should approach a constant different than zero at large times, in agreement with the DNS 
data. The RDT results suggest a persistent anisotropy of the higher odd order normalized transverse derivative 
moments at large Reynolds number, in agreement with the experimental findings^ 

There are numerous experimental investigations of incompressible homogeneous shear Q.ow^L^ however, no such 
study exists for the compressible case. Compressibility effects in turbulent flows are important in many practical 
applications ranging from combustion processes, to high speed aerodynamics, and to astrophysical phenomena. Al- 
though compressible turbulence has been the subject of intense research in the last 50 years, much less is known in 
comparison with incompressible turbulence^ In order to isolate the volume changes of the fluid elements, the velocity 
field is usually decomposed into a solenoidal (divergence free) part and a dilatational part. For homogeneous flows 
this decomposition is unique up to an additive constant, which can be taken to be zero without the loss of generality. 
Such decomposition in spectral space has been exploited by Rapid Distortion Theories in studies of shock-turbulence 
interaction and homogeneous turbulence subject to bulk compression or uniform mean shear ^1 

For the RDT equations of incompressible homogeneous shear flow, analytical solutions in wavenumber space are 
known for the velocity field in both viscous and inviscid cases42 The integration of these solutions over the wavenumber 
space can yield predictions for various quantities in real space. In general, only for very early and long times the 
integrals can be evaluated analytically. Long time solutions for the velocity variances and shear stress are obtained 
by Moffatt 13 and Rogers^ The early time evolution of the flow is correctly captured by the RDT equations since the 
non-normal amplification mechanism is linear in time. 15 ' 16 However, due to the neglecting of the nonlinear terms and 
thus the energy cascade to small scales, the RDT predictions can become eventually very different than those of the 
full nonlinear equations. Nevertheless, for the incompressible case, the RDT equations correctly capture the behavior 
of various correlations coefficients in the fully developed flow fields 

For the case of compressible homogeneous shear flow in the RDT limit no analytical solutions are known. Simone 
et al. 1 - performed RDT simulations of homogeneous shear flow and showed that the role of the distortion Mach 
number, M4, on the time variation of the turbulent kinetic energy is consistent with that found in the DNS results. 
They also identified different time regimes in which the various contributions to the terms in the RDT equations in 
spectral space change qualitatively, which might be the reason for the difference in the early and long time influence 
of Md on the kinetic energy growth. However, in the absence of analytical solutions, it is difficult to predict how each 
wavenumber will affect the result of the integration over the wavenumber space. For example, for the incompressible 
case, Rogers^ shows that at long times most of the contribution to the velocity variances and shear stress comes from 
a very narrow region in wavenumber space, which shrinks in time. 

The DNS results for compressible homogeneous shear flow of Blaisdell et al. — and Livescu et al. — reveal that the 
explicit dilatational effects tend to occur predominantly in the direction of the shear. Livescu and Madnia— examined 
the anisotropies of the solenoidal and dilatational motions for different values of the Reynolds and turbulent Mach 
numbers and found higher levels of anisotropy for the normal components of the dilatational dissipation tensor 
compared with those of the normal components of the solenoidal dissipation tensor. By examining the energy transfer 
leading to the anisotropy of the dilatational motions, Livescu et al. — showed that the nonlinear terms in the transport 
equation for the kinetic energy components do not have a significant contribution. This suggests that a linear 
mechanism might be responsible for the amplification of the dilatational motions in the direction of the shear. 

The present study thus aims: (i) to point out a lack of isotropy of the even order transverse derivative moments 
for both compressible and incompressible cases and show that this anisotropy is due to a linear mechanism, (ii) 
to examine the ability of the linearized equations to predict the persistent anisotropy of the normalized odd order 
transverse derivative moments found in the experimental data, and (iii) to clarify if the amplification of the dilatational 
kinetic energy and dilatational dissipation in the direction of the shear are also due to a linear mechanism. Moreover, 
the different levels of anisotropy of the small scale dilatational and solenoidal motions for the compressible case is 
also discussed. Additionally, analytical solutions are presented for the RDT equations of compressible homogeneous 
turbulent shear flow. 

The paper is organized as follows. Section II describes the governing equations and their linearization in the RDT 
limit. The numerical methodology for solving the full nonlinear and the RDT equations with the parameters for the 
cases considered is presented in section III. Section IV contains an analysis of the RDT equations for compressible 
homogeneous shear flow. The inviscid RDT equations are solved analytically in the incompressible limit {M to S — > 0) 
and pressure release limit (M to S — * 00). For finite values of M to S and large times, analytical solutions are presented 
for the Fourier modes of pressure and velocity components. In section V the analytical and numerical solutions of 
the RDT equations are compared with the DNS results. It is shown that the amplification of the dilatational kinetic 
energy in the direction of the shear and the anisotropy of the normal components of the dilatational dissipation 
tensor are captured by the RDT equations. Also in this section, the anisotropy of the even and odd order transverse 
derivative moments of the solenoidal and dilatational velocity fields for the compressible case and of the velocity field 
for the incompressible case are studied using DNS results and shown to be predicted by the RDT analysis. Summary 
and conclusions are given in section VI. 
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II. GOVERNING EQUATIONS 

The conservation equations which provide the mathematical model of the problem are the continuity, Navier-Stokes, 
and energy transport equations. For a compressible homogeneous shear flow (for which the mean, Favre averaged, 
velocity is given by U\ = Sx 2 with S constant, u 2 = 0, and M3 = 0) after Rogallo's transformation of coordinates^! 
x\ = BijXj, these equations, in their non-dimensional form, becom o 19 ! 22 

P,t + (pu'!), k B kj = (1) 

(pu"). t = - pu 2 SS l i - (pu"u"). k B kj - p, k B kl 

+ [(Tlj + ^iSSaS^ + SS^MkB^ (2) 

(p4>),t = S(r[ 2 - pu"u 2 ) - (pu"<f>), k B k j - (pu") !k B k j 

The stress tensor is t-j = ^^(s'^ — ^A'Sij), where s'^ — j{u'[ k B k j+u'- k B ki ) is the strain rate tensor, and A' = u" i k B ki 
is the dilatation of the velocity fluctuations. The primary transport variables are the density p, velocity fluctuations 
(with respect to Favre average) in xi direction u", and modified total energy <\> = j^~/n + \ u i u i \ where p is the 
instantaneous pressure and the ratio of the specific heats 7 = 1.4. The coordinate transformation matrix has constant 
diagonal components Bu = [3i and the only off-diagonal nonzero component is B12 = —(3\St. 

The thermodynamic variables are related through the equation of state, p — -^jti, and the nondimensional viscosity 

is modeled by assuming a power law p — T n , with n — 0.7. The reference scales used to nondimensionalize the 
governing equations are the initial rms velocity fluctuations (uq = \J u'/ u'/ ), initial mean temperature (T ), initial 
mean density (po) and a reference length-scale (Iq). Consequently, the non-dimensional parameters in Eqs. (0 are 
the computational Reynolds number, Rcq = Po ^°'° , the Prandtl number, Pr — and the reference Mach number, 

Mo = , where R is the gas constant and c p is the specific heat at constant pressure. The reference viscosity, 

Po, and thermal diffusivity, k , are assumed to be proportional to T ™. 

For the incompressible case, the energy equation is eliminated, the continuity equation is replaced by the incom- 
pressibility constraint, and the molecular viscosity is constant. The governing equations become 



l j,k 



B kj = (4) 



u i,t = -U2S6H - [u'ju'^^Bkj - p tk B kl + T' ij k B kj , (5) 
where u\ are the velocity fluctuations and = j^iu'i k B k j + u'j k B k i 
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A. Linearized equations for compressible homogeneous turbulent shear flow 

The linearized nondimensional equations of motion for density, p' , velocity, w", and pressure, p', fluctuations in a 
compressible homogeneous turbulent flow are derived by Kovasznay^ For the case of a homogeneous turbulent shear 
flow these equations correspond to the RDT limit, which assumes large values for the mean shear rate S. In the 
moving coordinate system defined above, the linearized equations become 



< p > A' (6) 

1 I 1-, r, lie 2 < p, > . . 1 

. —P,k B ki ~ Su 2 Su + - — — — (s - - -x 

< p > ' < p > Reo J 6 



u'lt = - ~ n ^ P[k B kt ~ Su 2 %-l + - : - p _ (s' i:j - -A 1 'Sij) ik B kj (7) 



, . , 4 < p > (7 - 1) „ , 

P'.t = -7<P>A' + ^— Ai >Ss' 12 

Re 

, ' /' V . (A- ) . (8) 



MyReoPr > \< P > < P> 
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While the average density is constant in time for homogeneous shear flow, the linearized equation for the mean 
pressure is 



<P>,= </X> ff (7 ~V (9) 
Re 

with the right hand side of the equation given by viscous dissipation of the mean flow. 

Although in general the linearized equations form a coupled system of equations, for the inviscid case they can be 
decoupled (note that in this case the mean pressure becomes constant). Thus, an equation for the dilatation can be 
obtained by taking the divergence of Eq. J7J) 



A> t = --L- p ' jkBjl B kt - 2Su'l k B kl . (10) 

This equation is used to eliminate the dilatation from Eq. © . Finally, after eliminating u'^ with the use of Eq. J7J , 
an equation for the pressure fluctuations is obtained 



^ = v\ m B ]% B kl - ASp{ 3k B 3l B k2l (11) 
c o 

where cq = (7 < p > / < p >) 1 / 2 (which for the inviscid case becomes constant). Equation (| 1 1 1) is a third order linear 
differential equation in time, with variable coefficients. For S — (decaying isotropic turbulence) Eq. (|1 l|l leads to 
a wave-like equation, with a bounded solution. For S 7^ 0, the time t can be nondimensionalized by 1/S, so that 
X] = St. Furthermore, after applying the Fourier transform, Eq. (|ll|l yields 



Aw? = ~^|( fc2 A'7 - 4fclfc2P), (12) 

where the wavenumber kj is related to the wavenumber fcj in the moving coordinate system by kj — B^ki, and 
k 2 = kjkj. The notations used in the present paper for the wavenumbers in the fixed and moving coordinate 
systems are the opposite of the notations used by Rogers^. Due to the time variation of k 2 , Eq. I|12(l does not 
have an easily derived analytical solution except in the case k\ = and for very early and long times (77 ~ 
and 77 — * 00). Nevertheless, it should be noted that, although for the full nonlinear equations the initial turbulent 
Mach number, M to — (2Kq) 1 ^ 2 /cq, and S are independent parameters, they appear together in Eq. 112|) . since 
cl/S 2 = 2K / < p > /(M to S) 2 = l/(M to S) 2 (note that the initial value of the turbulent kinetic energy is K = 0.5 
and < p >= 1). Moreover, if the distortion Mach number^ Md = SL/c, is introduced, then c^/S 2 — Lq/M^, where 

Md Q is the initial value of the distortion Mach number and Lq is the initial value of the integral scale L = (2K) 3 ^ 2 /e. 
Similar equations, depending on Cq/5 2 , rj, the wavenumber, and initial conditions can be found for the Fourier modes 
of dilatation, density fluctuations, and velocity components. In order to complete the analysis, the initial conditions 
for Eq. (|12fl are considered. These relations can be rearranged to isolate the parameter (M to S) by making the 
transformation p' — M ta p, 



p'\ v =o = M to p (13) 

p',r,U=a = -- — - ° Q do (14) 
' < p > M to b 

( k' 2 2k' \ 2K 

P '^=° = {-jMZW {MtJo) + Tlts iiL2 ° ) 77^' (15) 

where k\ is the initial value of the wave number in the fixed coordinate system (k\ — hi), k' 2 — k'^k'^ i — \f— 1, 
and do and u 2a are the initial values of the Fourier modes of the dilatation and velocity fluctuations in x 2 direction, 
respectively. It can be seen that, if the initial value of the pressure mode scales with l/M to for constant (M to S), then 
p' depends on (M to S) and rj only. This condition is fulfilled if the pressure fluctuations are initialized with zero or 
from a Poisson equation, which is the case with most DNS studies of homogeneous turbulent shear flow. Note that 
to first order the Poisson equation becomes 



M to p = 2i <p> (M ta S)-±u 2 . 



(16) 
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In this case, from the inviscid momentum equations and Eg. 1101 for the dilatation it yields, after replacing t with 77, 
that Hi and A depend only on (M to S). Moreover, if the initial values of the density fluctuations are set to zero, Sp 
will also be a function of (M to S), wavenumber, and rj only. 

As a result, for the initializations usually used in the DNS of homogeneous turbulent shear flow, the inviscid 
linearized theory predicts the existence of only one compressibility parameter, (M to S), for the properly scaled variables. 
However, the addition of either the nonlinear or the viscous terms prevents the governing equations to be written 
in terms of only one compressibility parameter. Therefore, for moderate values of S or Reynolds number, two 
independent compressibility parameters can be defined. Nevertheless, as pointed out by Simone et al. j 17 changes in 
Md at constant M to are more important to the growth of turbulent kinetic energy than changes in M to at constant 
M do . 

III. NUMERICAL SOLUTION PROCEDURE 

A. DNS methodology 

Equations 0-0 are solved using the Fourier pseudo-spectral method with Rogallo's remeshing technique,— Details 
can be found in Livescu et al. i 9 All simulations are performed within a box containing 256 3 points. The computational 
domain is twice as long in the stream-wise direction as in the cross-stream and span-wise directions so that fj\ — 0.5, 
/?2 = 1-0 and /?3 ~ 1.0. The initial velocity is a random solenoidal field with unity r.m.s. and the spectrum given by 
E{k) = fc 4 //&, exp(-2k 2 /k 2 v ) with k 0v = 10. The initial pressure fluctuations are evaluated from a Poisson equation, 
the initial density field is uniform, and the initial value of the mean pressure is calculated from the mean equation 
of state. The numerical method, mesh size and initial conditions for pressure and velocity fields are the same for the 
incompressible case. 

Several cases, 1 to 5, are considered at different values of the mean shear rate. For each case two simulations, 
one compressible and one incompressible, are performed. Table I presents the relevant information for each of the 
cases considered. The initial nondimensional shear rate is defined as Sq — 2KS/e. For the compressible cases, 
the distortion Mach number Md is also listed. Cases 1 to 5 have the initial value of the Taylor Reynolds number 
Re\ = 21 and it increases to Re\ ~ 100 at the end of the simulations. Additionally, one more case (case 6) is 
considered with Re\ = 50, and all the other parameters as those of case 2. In the next sections, the compressible 
simulations corresponding to cases 1 to 6 are labeled CI to C6, and the incompressible cases are labeled II to 16. 
All simulations were monitored to ensure that the integral scales remain small compared to the box size and the 
Kolmogorov microscale is larger than the grid size. The compressible runs were also repeated with initial turbulent 
Mach numbers ranging from M to — 0.1 to M ta — 0.6 and all the results presented in the next sections remained 
qualitatively unchanged. 

Based on the magnitude of the shear rate, a simple scaling argument can be used to characterize the relative 
importance of the nonlinear term in the momentum equations compared to the linear term arising from the presence of 
the nonzero mean velocity gradient. Thus, as shown in Ref. l24L the nonlinear term is negligible when S*/Re\ 1 3> 1/8 

and dominates when S*/Re\ 1 -C 0.35, where S* — u' 2 S/e and Re\ x = u' 2 /(v^J < u\ l >), with u' = ^<u\>. 

Figure 1 shows the time evolution of Sl/Re\ 1 for the compressible cases considered. The results obtained for the 
incompressible cases are close. The initial range of variation covers the region in which both the linear and nonlinear 
terms in the momentum equations are important, with case CI, with lowest shear, for which the nonlinear effects 
dominate and case C5, with highest shear, for which the linear term becomes more important. After a development 
time the results obtained for the cases with different values of the mean shear rate become close, while retaining the 
dependence on the initial Reynolds number. For all cases, the late time values of Sl/Re\ 1 are in the region where 
the nonlinear effects dominate, based on the scaling arguments given in Ref. [3 

B. RDT methodology 

The unknowns in Eqs. 0-101 are transformed intoii 
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fi(k,r]) = 17) 

TO 

fa(k,T]) = (klUl + k 3 U 3 )—?- - ^U 2 (18) 

; \ fciui + k 2 u 2 + k 3 u 3 A 

MM) = £ = ^ (19) 

/ 4 (M) = »P (20) 

/ 8 (M) = V> (21) 



where to = \/k\ + fef. /i is the component of the velocity along the perpendicular on the plane defined by the 
wavenumber vector and fc 2 direction, fa is the component along k 2 direction so that /i and fa form the solenoidal 
velocity and fa is the projection of the velocity on the wavenumber vector and corresponds to the dilatational part 
of the velocity. Additionally, since the equations arc linear, a matrix exponentiation method can be used by letting 
fi(k,rj) = g ij (k,ri)f j (k,0) ) with g tj (k,0) = 5 y . 

The resulting system for the unknowns gij becomes 

93j ~ p „ ~ — Q k 9i 3 (22) 





k 3 I 


9ij,v 


= k 92j + ~ 




kik 2 


92j,V 


k 2 92, ■ 




k\m 


93j,V 


2 k 2 



mk Reo < p > S 

k 2 g 2j (23) 



ki _ < p > , 2 



Ti^ J Re < p> S 
kik 2 k 4 <p> 2 

T72"^3j - - — — ^ - o~5 — ~ — TT; fc 9zj (24) 
kk z < p> b 6 Re < p > b 



7 < p > 2<p> ( kiks k\k 2 n k\k 



9«,n - ~^ k 9^ + -r^-(7 - 1) + feJTflW + 

< P > , 2 ( 94j 95j 



M?Re PrS \<p> < p > 



(25) 



g5j,r, = < ^, > fcg3j, (26) 

where j = 1,5. The svstem F22l26l is transformed into spherical coordinates by setting k[ — k' cos a, k 2 = k' sin a sin (f>, 
and fc 3 = fc'sinacos^. Cases with (101 x 602 x 204) grid points for the coordinates (k',a,4>) are considered. The 
discretization is nonuniform, allowing for more grid points near the line a = <f> = tt/2. This is necessary, as explained 
in the next section, in order to correctly capture a local peak of some of the unknowns, which gives an important 
contribution to their integral over the wavenumber space. A fourth order Runge-Kutta scheme is employed for the 
time advancement. The equations are also solved for the incompressible case (which is recovered by letting fa — 0). 
It should be noted that the system 1221261 is independent of the initial spectra, so that it is solved only once for a 
particular set of parameters. Furthermore, when integrating over the wavenumber space only the initial power spectra 
are required, and not the phase information. 

The linearized Eqs. (|22|I - H26|I are solved for different values of Sq ranging from 7.24 to 362. For all compressible 
cases M to = 0.3, so that the range in Sq corresponds to 2.2 < Md < 109. This range extends from a low shear 
case in which the nonlinear effects in the corresponding DNS case are important to cases with very high shear. The 
initial velocity power spectrum is the same as in the DNS simulations described earlier. For the compressible case, 
for each value of (M to S) four types of initial conditions, labeled LI to L4, are considered (Table II). For cases LI and 
L3 the pressure fluctuations are initialized from a Poisson equation, while for cases with L2 initial conditions they 
are set to zero. The initial velocity field is solenoidal for cases LI and L2 while xo — Kdo/Ko = 0.1 for cases L3 
and L4, where Kdo is the initial dilatational kinetic energy. For cases with L4 initial conditions the strong form of 
the acoustic equilibrium is considered^ so that there is an equipartition of energy between the potential and kinetic 
energies of the dilatational component at each wavenumber. Therefore the initial pressure power spectrum is given 
by E p {k') = y/j < p > | t= o < po > Ed(k'), where Ed is the initial power spectrum of the dilatational velocity. 
For each case, both the inviscid and viscous equations are solved. The types of initial conditions considered are not 
exhaustive for the initial conditions used in DNS of compressible turbulence (e.g. see Refsi^2& for decaying isotropic 
turbulence). However, they were chosen such that the role of a nonzero initial dilatational field can be isolated. In 
the next sections, the numerical solutions of the RDT equations are compared with the DNS results and with the 
analytical predictions. 
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IV. ANALYTICAL RESULTS 

By examining Eq. l(T2l for the pressure mode three cases can be identified based on the values of l/{M ta S) 2 . 

A. "Pressure- released" limit ( , , . 1 „..> — * 0) 

This case corresponds to values of the distortion Mach number approaching infinity. After integrating Eq. I|12[) . a 
quadratic dependence on r\ is obtained. However, when the initial conditions (|13|) - (|15[1 simplified for large (M to S) are 
applied, the solution is obtained as p = p\ v =o. Since the pressure mode is multiplied by the factor l/(M to S) in the 
transport equations for the velocity modes, it can be dropped from these equations. Therefore the equations for the 
velocity modes become: 

Ml, 77 = ~U2, U 2 ,r, = 0, U 3:V = (27) 

with the solutions u\ — U2|?7=o?7 + u\\ n= a, ii2 — "2^=0, and U3 = M3^ = o- After integrating over the wavenumber 
space, it is obtained that the kinetic energy in X\ direction increases quadratically, while it remains constant in the 
other two directions. This case corresponds to the "pressure-released" limit discussed by Cambon et al. — and Simone 
et al. il. 

B. Incompressible case ( pg~ s yi — * 00 ) 

Since S is assumed large in the RDT limit, as (jv/ t 1 g)2 — > 00 it yields M to — ^ 0, so that the incompressible case is 
recovered. The transport equation for the pressure mode becomes 

k 2 p^ - 4fcifc 2 p = 0. (28) 
The solution is straightforward and can be written as 

k' 4 

P = p\ v =o-^, (29) 

with p\r]=o = c 2iS-piU2\-q=Q- In order to find the pressure variance from Eq. I|29() . the integration over the wavenumber 
space should be performed numerically Analytical expressions can be obtained to approximate the time behavior of 
the pressure variance only for 77 — > (early time) and r\ — > 00 (long time). In both cases the integration is conveniently 
done in spherical coordinates. The procedure is similar to that used by TownsencU^ and Rogers^ to find the time 
dependence of velocity and scalar one point statistics in incompressible turbulent shear flow. After the transformation 
k[ = fc'cosa, k' 2 — k' sin a sin </>, fc 3 = k' sin a cos </>, and dk' = fc' 2 sina da d<fi dk' , for isotropic initial conditions (so 

that I "0.2 1 77=0 1 2 = rri fcl fc !if 3 with E(k') the initial velocity power spectrum), the pressure spectrum function becomes 



E(k') AS 2 cos 2 a(cos 2 a + sin 2 a cos 2 * 
4irk' 2 k' 2 (l — 2rj cos a sin a sin <f> + rj 2 cos 2 a) 4 ' 



. / . j uud ut^uo uc ~r 0111 utua yj f , , 



and the pressure variance can be found from 

<p ,2 >(v)=^- r * P dk[. (31) 

For small values of rj (early time) the integrand can be expanded in a Taylor series in 77, which can be integrated 
analytically. It yields 



<p' 2 > 4 , , 



. — V +-, (32) 
15 105 ' y ' 
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where q — AS 2 J E ^ dk'. In terms of the initial value of the pressure variance the dependence is ^yg^ = 1 + |?? 2 + ... 

As explained by Rogers^* the long time behavior can not be obtained using the same procedure, since expanding 
around 77 = oo yields a result which is not uniformly valid over the entire domain. This can be easily seen by evaluating 
p for k' 2 = 1, &3 = 0, and k' x = I/77, which for 77 — > oo yields 



lim p(k[ — I/77, 77) 



(33) 



Therefore a significant contribution to the integral of the pressure spectrum over the wavenumber space comes from 
a narrow region near k[ = which decreases in size as 77 increases. Following Rogers; 1 - the asymptotic behavior can 
be found by dividing the domain into two regions, an inner region located around rjk^ = 0(1) (which corresponds 
to a = 7r/2) and an outer region consisting of the rest of the domain covered by < a < ir and < <fi < 2tt. In 
the inner region, around a = tt/2 and <f> = 7r/2, $ p (fe')fc' 2 sina has a peak. The magnitude of the peak and the 
extent in a and <fi over which it occurs can be found by writing a = tt/2 — e a and 4> — tt/2 — and expanding the 
sin and cos functions in the expression for the pressure spectrum. It yields that <& p (fc')fc' 2 sin a is of order rf for 
e a = I/77 ± 0{l/r} 2 ) and = 0(1/77), so that the peak occurs over a region of extent 0(l/n 2 ) in a and 0(1/77) in 
fa This indicates a linear increase for the pressure variance. The correct behavior can be obtained by combining the 
results from the inner and outer regions. The result obtained remains < p' 2 > / < p$ >~ 77, also supported by the 
results of the numerical integration. 

C. ... 1 .ii finite 

In this case the transport equation for the pressure mode (Eq. I12|) remains third order, and does not have an easily 
derived analytical solution. However, some limiting cases can be found. First, it should be observed that p approaches 
zero as k[ — > 00 so that only finite values of k[ are of interest. Then for small values of 77 the equation reduces to an 
equation with constant coefficients which can be easily solved. However, since the roots of the characteristic polynomial 
are dependent on the parameter (M to S), the analytical computations involving the integration of the solution over 
the wavenumber space are cumbersome. Nevertheless, the solution was verified numerically and it agrees reasonably 
well with the full numerical solution of the RDT equations. 

As 77 — + 00, unlike in the incompressible case, the behavior of the Fourier modes in the outer region (k^rj >> 1) 
is very important for the integration over the wavenumber space. It is shown above that for the incompressible case 
the contribution from the inner region yields a linear increase for the pressure variance. However, as the parameter 
(Ah s)^ decreases, the peak of the pressure mode occurring in the inner region decreases its magnitude (Fig. 2) and 
the contribution from this region to the pressure variance yields a less than linear increase in time. This contribution 
to the pressure variance is negligible at long times since it is shown below that the contribution from the outer region 
leads to a linear increase in the pressure variance. 

In the outer region ^77 — > 00 for 77 — > 00 and Eq. I|12f> becomes 



2K k' 2 



(v 2 P,v +4t7p). 



(34) 



< p > {M t0 S) 2 

The solution of Eq. 134(1 can be written as p = cf>i p\ri=Q + 4>2 P,?i|?7=o + 03 P,r/ri\ri=o, with fa given by the formulas 



1F2 [1] 



1 3 

2'4 



4 s 



.(C)(C/2) 3/4 2r( i 



M t0 S 



1/2 



< p > 

2K 



<P> 
2K 

1/2 



1/4 



(35) 
(36) 
(37) 



where 1F2 and J v represent a generalized hypergeometric function and the Bessel function of the first kind of order 

k[v 2 I K a 



v 1 respectively, T is the gamma function and £ 



M t0 S\ 2 <£> 

the asymptotic expansions for these functions given in Ref. |23, it is obtained that 



For large values of the nondimensional time 77, using 
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_ (C/2) i/4 sm(C _^ )r (-) (38) 



V2 „ ^ \ 1/4 



fc „ (c/2) v^ 8(c _ i)2r (|)(^) (39) 
*»«« 1/4 -«-f) r (l)^|(^) 1/2 - « 

These expressions indicate that the Fourier mode of the pressure oscillates in time and the amplitude of the oscillations 
increases as ^/rj. In order to find the pressure variance, the integration over the wavenumber space can be performed 

using the Riemann-Lebesgue theorem, f(k')e lk n 6.k' — > as ?y — * oo provided that J \f(k')\dk' exists. Since 
sin 2 a = (1 — cos2a)/2 and cos 1 a = (1 + cos2a)/2, the integration yields a linear increase in time for the pressure 
variance. This is consistent with the numerical solutions of the linearized inviscid equations presented in Fig. 3. 

Similar to the pressure variance, it can be shown that for finite values of jjj^-gjz, most of the contribution to the 
dilatation variance and kinetic energy in xi direction, K2, comes from the outer region. From the solution obtained 
for the pressure mode, approximations for the asymptotic behavior of the dilatation and 112 velocity modes in the 
outer region can be derived using equations 



P'r, (41) 

2K ' ' 

- = i\v ik 2 f>' , , 

2Ai 2 < p > (M t0 S)h' 



The relations obtained are A = (j> 1A p'\ v =a + <p2 A f, v \n=o + h A p[ m \v=v and "2 = 4>iv3, f\v=o + 4>iui p[ r ,\v=o 
hu2 p[r,ri\v=o witn ^iA and given by 




Again, using the Riemann-Lebesgue theorem, Eqs. (|43[) - l|48[) can be integrated over the wavenumber space to yield 
an rf increase in time for the dilatation variance, < A' 2 >, and a linear increase for < u\ >■ This represents an 
important difference compared to the incompressible case, where it it is known that for large times the u-i velocity 
variance decreases in time as In(4?7)/(4?7)J4 It is shown below that the increase in time in the compressible case is 
due to the dilatational component. 

In the incompressible case, the influence from the inner region to < u\ > and < it§ > is significant and at large times 
< u\ >« 2 In 2 77 and < u\ > ir 2 /8\ni] — Cw^ For the compressible case, after using the approximation ki ~ —tjk^, 
the transport equations for u\ and U3 can be solved analytically and for large values of 77 
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mi ~ u 1 \ v=0 + — sin(—— g\hr- — - + d2 (49) 

y/n M to s y 2 < p > 



u 3 ~ «sU=o + ^ s -(M^\/2^ + e2) (50) 

with the coefficients <ij and e< depending on (M to S) and the wavenumber. Therefore, in the outer region, for large 
times, iii an d u 3 become constant. Numerical solutions of the RDT equations indicate that the contributions from 
the inner region to < u\ > and < u 3 > remain dominant for the compressible case and K\ and K 3 increase in 
time, as shown in Fig. 4. However, as the compressibility parameter M to S increases, these contributions decrease 
their importance resulting in a reduced rate of increase of the turbulent kinetic energy in x\ and X3 directions. This 
explains the "stabilizing" effect of compressibility on the evolution of the turbulent kinetic energy found in the RDT 
study of Simone et al. H 

From the solution corresponding to the Fourier mode of the dilatation, the dilatational parts of the velocity modes 

can be obtained as iii d = — t^A. Numerical solutions of the RDT equations indicate that the inner region has a 
negligible contribution to the variances of the dilatational velocity components. In the outer region, for large values 
of 77, the time variation of the dilatational parts of the velocity modes can be derived from formulas 1431451 as 



en . k[rf / K 
^ ~ ^ Sm( M^\/2^ + Cl2) (51) 

U2 d ~ c 2 iy/rj sin(——*l "" + c 22 ) (52) 




y/rj M to b V 2 < p > 



where the constants depend on (M t(l S) and the wavenumber. 



V. COMPARISON WITH DNS 



Although the early time response of the flow to the presence of the shear is correctly captured by the linearized 
equations, the long time evolution of the flow can become very different than that obtained from the full nonlinear 
equations. However, there are quantities which follow the RDT predictions even at long times. This is the case with 
various correlation coefficients (e.g. the Reynolds stress correlation coefficient) in the incompressible case^ For the 
compressible case, the effect of compressibility on the evolution of the turbulent kinetic energy is correctly captured 
by the linearized equations^ 

It is shown below that other characteristics of the dilatational field are captured by the linearized equations. In 
particular, it is shown that the amplification of the dilatational kinetic energy in the direction of the shear and the 
anisotropy of the normal components of the dilatational dissipation rate tensor are captured by the RDT equations. 
Furthermore, DNS results show that the transverse even order derivative moments of the velocity field in both 
compressible and incompressible cases are anisotropic and the anisotropy increases with the order of the moment. 
Thus, the anisotropy of the higher order moments may persist at large values of the Reynolds number, in contrast 
with the local isotropy principle, — Moreover, the DNS results indicate that the odd order normalized X2 derivative 
moments of u± approach a constant different than zero at large times. These results can also be explained using the 
RDT solutions. 



A. Dilatational kinetic energy and dissipation anisotropy 

By comparing formulas 15 11531 with solutions obtained for the velocity modes, it can be seen that at long times 
iii ~ Ui a , iii ~ U2 d , and ii 3 ~ u 3a . It is obtained that most of the dilatational kinetic energy resides in the direction 
of the mean velocity gradient, consistent with the numerical solution of the linearized Eqs. 10-1© (Fig- 5a). The 
addition of the viscous terms decreases the dominance of K2 d component compared to the other two components, 
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as can be seen in Fig. 5(b). A further decrease of the relative magnitude of K% d is obtained in the fully nonlinear 
case (Fig. 5c). Nevertheless, as Fig. 5(c) shows, the dilatational kinetic energy in the direction of the mean velocity 
gradient obtained in DNS is larger than in the other two directions, in agreement with the previous studies £LliL2iL 
The above analysis indicates that the amplification of the dilatational kinetic energy in the direction of the mean 
shear can be explained by linear effects. 

Previous numerical simulations of compressible homogeneous shear fiowiSii^iSll indicate that the dilatational dissi- 
pation rate tensor behaves very differently than its solenoidal counterpart. Thus, the dilatational dissipation rate is 
amplified in the direction of the mean shear. This behavior was previously attributed to the formation of shocklets 
aligned preferentially in X2 direction. Livescu and Madnia 20 performed simulations at different values of the turbulent 
Mach number and obtained no significant decrease in the anisotropy of the dilatational dissipation rate tensor in low 
turbulent Mach number simulations. Therefore, the amplification of the dilatational dissipation rate in the direction 
of the mean shear is not associated with the presence of shocklets. This conclusion is also obtained in Ref. 0. Figure 
6 shows that 622 d becomes more important compared to en d and €33 d as M ta S, which is proportional to the distortion 
Mach number, increases. After multiplying the solutions for the dilatational velocity modes (Eqs. 151153(1 by iki and 
performing the integration over the wavenumber space, it can be seen that in the RDT limit < Q^t ^ becomes 

much larger than < -g^ > and < ^f d >. Numerical solutions of the RDT equations in both inviscid and 
viscous cases confirm this result. Therefore, the anisotropy of the normal components of the dilatational dissipation 
rate tensor can be associated with a linear mechanism. 



B. Higher order derivative moments 

In the outer region &2 increases continuously, while it is very small in the inner region (where kl Y r\ ~ k' 2 ). Using the 
solutions provided in the previous section for the solenoidal and dilatational components of the velocities, it can be 
shown that after the multiplication by &2, the contribution from the outer region becomes dominant at long times. 
Therefore, the X2 derivatives of the velocity components become much larger than the derivatives in the other two 
directions. As a result, all the even order derivative statistics involving X2 derivatives become amplified and this effect 

is more pronounced for higher order statistics. Figure 7 compares moments of the transverse derivatives, (i -^f: J j 
where j = 2,3 and n — 2, 4 and 6, of the dilatational velocity in x\ direction. It can be seen that, as the distortion 
Mach number increases and the RDT limit is approached, the moments containing X2 derivatives become much larger 
than the corresponding x% derivative statistics. Moreover, in agreement with the RDT predictions, this anisotropy is 
strongly amplified for higher order moments. This again suggests that the small scale anisotropy of the dilatational 
motions found in the DNS results is produced by a linear mechanism. On the other hand, the results presented show 
that the higher order transverse moments are increasingly more anisotropic as the order of the moment increases. 
Thus, even if the anisotropy decreases with Reynolds number, some degree of anisotropy will persist for higher order 
transverse derivative moments at large values of the Reynolds number. 

Although the inner region has a significant contribution to the one-point statistics of the solenoidal velocity com- 
ponents in x\ and X3 directions, this contribution becomes negligible for the statistics involving X2 derivatives. As 
explained above, in the outer region u\ g and u^ s become constant at large times so that after the multiplication 
with k2 they increase linearly. Therefore, the outer region yields a quadratic increase in time for < g" 1 ^ > and 
< d"3 3 du 3s ^ rp^ c con tribution from the inner region does not increase after the multiplication with fc 2 and the 
inner region yields a less than quadratic increase for the above quantities. Similar results are obtained for higher 
order moments. Thus, similar to the dilatational velocities, in the RDT limit the even order transverse derivatives of 
the solenoidal velocities in x± and X3 directions are anisotropic and become even more anisotropic as the order of the 
statistics increases. Figure 8 confirms that this effect is also present in the DNS results for the derivative moments of 
Ui s and similar results are obtained for the derivative moments of U3,, . 

A comparison between Figs. 7 and 8 indicates that the transverse solenoidal derivative moments are less anisotropic 
than their dilatational counterparts. Again, it can be shown that this is also in agreement with the RDT predictions. 



Thus, since Ui B becomes constant in the outer region, it yields that k2U\ s ~ rj and / J ~ V 2 - However, the X3 

derivative retains influence mostly from the inner region and / ( ° ^ \ ~ 17. On the other hand U\ d retains influence 



mostly from the outer region. The Riemann-Lebesgue theorem can be used to integrate formula 1511 multiplied by 
and fc 3 to obtain for the dilatational component that ^ ^ ^dxt ) j ~ ^ an< ^ ( ~5xi ) ^ ~ ^/v- Similar results are 
obtained for higher order statistics. Therefore, in the RDT limit, the even order transverse derivatives of the u\d 
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are more anisotropic than their solenoidal counterparts for large but finite values of rj (Fig. 9). Similar results are 
obtained for the velocity in X3 direction. 

Using the same analysis as above, it can be shown that the odd order moments of the xi derivatives of u\ s also 
have most of the contributions from the outer regions. Thus, after the multiplication with fcj ~ ( — V^'if 1 -! the solution 
in the outer region becomes much larger than that obtained in the inner region as r\ — > 00. Note that \u\ s \ increases 
in time as rj 1 / 2 in the inner region and the multiplication by &2 does not increase this rate. Since in the outer region 
ui a becomes constant at large times, it yields that the moment of order n increases in time as rf 1 , and, therefore, 



the normalized odd order moments, 



with n odd, should approach a constant different 



than zero at large times. As Fig. 10 shows, this prediction is in general agreement with the DNS data. The same 
analysis can not be used for u^ 3 . It was shown that in the outer region v,3 s approaches its initial value. Therefore, 
after multiplication by k% ~ (— Tjk'i)'" the odd order moments of the x-i derivative, would approach values proportional 
to those of the initial odd order moments of the x\ derivative. However, for isotropic initial conditions, these values 
are zero. This is consistent with the overall symmetry of the problem, which requires that the odd order x\ and xi 
moments of the U3 velocity remain zero at all times. 

For the incompressible case, complete analytical solutions are known for the velocity fieldiiSiii Although most of 
the contributions to the velocity variances in x\ and X3 directions come from the inner region, it is easy to show 
that the statistics involving X2 derivatives have contributions mostly from the outer region, where at large times 
fc 2 ~ — Tjk[. Moreover, the behavior of the velocity modes in the outer region is similar to that obtained above for 
the solenoidal velocity modes. Consequently, in the RDT limit the transverse even order moments of the velocity 
derivatives are anisotropic, and the anisotropy increases for higher order moments. Similar to the solenoidal moments 
for the compressible case and consistent with the RDT predictions, the incompressible DNS results presented in Fig. 
11 show an increase in the anisotropy for higher order transverse derivatives moments of u\. Similar results are 
obtained for U3. Therefore, it is expected that for higher order moments the anisotropy will persist at large Reynolds 
numbers. 

In addition, similar to the solenoidal velocity for the compressible case, the RDT solutions for the incompressible 
case predict that the odd order normalized X2 derivative moments of u\ and 113 should approach a constant different 
than zero at large times. The DNS results are similar to those shown in Fig. 10 for the solenoidal velocity field in 
the compressible case and support the RDT prediction. As explained above, as the Reynolds number increases the 
small scales become more energetic and the flow departs from the RDT limit (see also Fig. 1). For very large values 
of the Reynolds number, the postulate of local isotropy requires that the odd order transverse derivative moments 
should approach zero. However, since the linear effects lead to a rf 1 increase in time for the xi derivative moments, 
for higher order moments the linear contribution may become dominant at long times even at large values of the 
Reynolds number. Therefore, the higher order normalized transverse X2 derivative moments could have persistent 
non-zero values at large Reynolds numbers. The present numerical simulations do not have a large enough Reynolds 
number to examine this hypothesis, however the analysis presented is consistent with the experimental results of Shen 
and Warhafti At high Reynolds numbers, the experimental results indicate that lower odd order moments of 
become small as the Reynolds number increases. However, this tendency is not observed for higher order moments. 
As suggested by Schumacher et al. ^4 the small-scale intermittency might play a role on the observed behavior of the 
higher order moments. Nevertheless, the persistent anisotropy of the higher odd order normalized derivative moments 
found experimentally is in agreement with the present linear analysis. 



VI. SUMMARY AND CONCLUSIONS 



The structure of homogeneous turbulent shear flow is studied using data generated by DNS and an RDT analysis 
for both compressible and incompressible cases. For the compressible case, simulations with different initial values 
of the distortion Mach number, M4, and Reynolds number are considered. Incompressible simulations are performed 
with the same initial pressure and velocity fields as for the compressible simulations. 

Previous DNS studiesi&iiLSi indicate that there is an anisotropy among the dilatational kinetic energy components. 
Moreover, the results of Livescu et al. — suggest that the energy transfer through the nonlinear terms in the transport 
equations for K% d is not responsible for this behavior. In order to verify that the anisotropy of the dilatational kinetic 
energy components can be explained by linear effects, the linearized equations are considered. The RDT equations 
are solved numerically, for both the inviscid and viscous cases, for a large range of distortion Mach numbers and 
different types of initial conditions. The RDT results are consistent with the DNS findings for the development of the 
dilatational field. 

Furthermore, for large times, analytical solutions are found for the inviscid linearized equations in Fourier space. 
The integration of these solutions over the wavenumber space can predict the behavior of various statistics in real 
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space. It is shown that the solutions for the pressure, dilatation and dilatational velocities Fourier modes for finite 
values of k[ oscillate in time and the amplitudes of oscillations increase. The analytical relations indicate that the 
magnitude of the Fourier mode of the dilatational velocity in the direction of the mean shear becomes much larger 
than in the other two directions. 

It is known that for incompressible homogeneous turbulent shear flow in the RDT limit most of the contributions 
to the velocity variances come from a narrow region in the wavenumber space situated near k[ = 0^ A similar trend 
is found for the pressure variance. However, it is shown that the contribution from this region to the statistics of the 
velocity field involving x-i derivatives for both compressible and incompressible cases and to all the statistics of the 
dilatational velocities becomes negligible at long times. For finite values of k[, the wave number in the direction of 
the mean shear increases continuously in time. Therefore, the derivatives in this direction become much larger than 
in the other two directions. As a result, in the RDT limit the even order transverse derivative moments of both the 
dilatational and solenoidal velocity fields for the compressible case and the even order transverse derivative moments 
of the velocity field for the incompressible case arc anisotropic. Moreover, the anisotropy becomes amplified for higher 
order moments. The RDT analytical solutions also predict that the small dilatational scales are more anisotropic than 
their solenoidal counterparts, in agreement with the DNS results. Since the X2 derivatives of the velocity field retain 
influence mostly from the region where the wavenumber in x% direction increases in time, it is shown that in the RDT 
limit the odd order normalized X2 derivative moments of u\ for the incompressible case and u\ s for the compressible 
case should approach a constant different than zero at large times. This prediction is consistent with the DNS results. 

As the Reynolds number increases the small scales become more energetic and for constant mean distortion the flow 
departs the RDT limit. Therefore, it is expected that, consistent with the postulate of local isotropy, the anisotropy 
of the small scales should decrease at large values of the Reynolds number, i.e. the odd order normalized x% derivative 
moments approach zero and the even order transverse moments become isotropic. However, it is argued that, since 
in the RDT limit the odd order moments increase in time as 77™, where n is the order of the moment, for higher order 
moments the linear contribution may be important even at large values of the Reynolds number. Thus, the RDT 
analysis suggests a persistent anisotropy of the higher odd order normalized derivative moments. In addition, since 
the anisotropy of the even order transverse derivative moments increases with the order of the moment, it is expected 
that persistent anisotropy levels in the transverse higher even order derivative moments will be found at large values 
of the Reynolds number. Although high Reynolds number results will not be available from DNS in the near future, 
they are within the reach of experiments and these predictions can be verified using experimental data. 
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Case 


Re\ 




Afj 


1 


21 


4.86 


1.46 


2 


21 


7.29 


2.19 


3 


21 


9.71 


2.91 


4 


21 


14.6 


4.37 


5 


21 


21.9 


6.56 


6 


50 


17.4 


5.22 



TABLE I: Parameters for the DNS cases. ' Values corresponding to the compressible cases. All compressible cases have 
Mto = 0.3. 
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| Type of initial conditions 


P\v=o 


Xo | 


LI 


Poisson equation 





L2 








L3 


Poisson equation 


0.1 


L4 


Acoustic equilibrium 


0.1 



TABLE II: Initial conditions for the linearized equations. 
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Figure captions 



Fig. 1 Time evolution of S* /Re\ 1 . 

Fig. 2. Compressibility effects on the k' x variation of the pressure mode. Inviscid case, r\ = 10, &2 — 1> ^3 = 0, 
with initial conditions p|^ = o = 1, P.rj\ri=o = 0, and p. nv \ v =o — 0. (a) Incompressible, (b) compressible, M ta S — 0.01, 
and (c) compressible, M to S = 1. 

Fig. 3. Numerical solution of the inviscid RDT equations for the pressure variance for M ta S ~ 4.6 (thin lines) and 
M to S = 31 (thick lines) and different initial conditions. For clarity the results obtained for M to S = 4.6 are magnified 
5 times. 

Fig. 4. Compressibility effect on the kinetic energy in x\ and £3 directions. Inviscid case, with LI initial conditions. 

Fig. 5. Time evolution of the normalized dilatational kinetic energy components for M to S — 1.53, 3.1, and 4.6. 
(a) Inviscid RDT, (b) viscous RDT, and (c) DNS. The curves correspond to cases C2, C4, and C5. 

Fig. 6. Time variation of the diagonal components of the dilatational dissipation rate tensor for M to S = 1.53, 3.1, 
and 4.6, corresponding to cases C2, C4, and C5. 

Fig. 7. Even order transverse derivative moments of the dilatational velocity in X\ direction, M% y = 



and Mt = ... XV 3 /,/ vn ■ (a) n=2, (b) n=4, and (c) n=6. 
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Fig. 8. Even order transverse derivative moments of the solenoidal velocity in x\ direction, M^ y 
■ flm \ "t" 2 / 7 i3ui.\n\ ■ n=2 (no symbols), n=4 (open symbols), and n=6 (closed symbols). 



Fig. 9. RDT results for the second order transverse derivatives of u\ s and u\ d . LI initial conditions, M to S = 1.53. 
Fig. 10. Odd order normalized transverse derivative moments of the solenoidal velocity in x\ direction, S ny = 

(fen 



. a) n=3 b) n=5. 



(((*)"»' 

Fig. 11. Even order transverse derivative moments of the velocity in x\ direction, incompressible case. M ny 
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